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A new method of cloaking is presented. For two-dimensional quasistatics it is proven how a single active 
exterior cloaking device can be used to shield an object from surrounding fields, yet produce very small scattered 
fields. The problem is reduced to finding a polynomial which is approximately one within one disk and zero 
within a second disk, and such a polynomial is constructed. For the two-dimensional Helmholtz equation, it 
is numerically shown that three active exterior devices placed around the object suffice to produce very good 
cloaking. 

PACS numbers: 42.25.Bs, 43.20.+g, 41.20.Cv 



Making a body truly invisible, in the sense of preventing 
it from absorbing radiation or scattering radiation in any di- 
rection, has long been regarded as the domain of science fic- 
tion rather than science, with a few exceptions |Q] |2|. This 
perspective has changed in recent years due to fascinating ad- 
vances in our understanding of how electromagnetic fields can 
be manipulated. Recent proposals for invisibility, as reviewed 
in lOJIU, can be broadly placed in two main groups: interior 
cloaking, where the cloaking device surrounds the object to 
be cloaked, and exterior cloaking, where the cloaking region 
is, surprisingly, outside the cloaking device. 

Interior cloaking methods include plasmonic cloaking due 
to scattering cancellation [5 |, transformation based cloaking 
16] El El HI, and active cloaking ifTUl . Transformation based 
cloaking has gained particular attention, being supported by 
experiment ifTTI fl2l [PJl [141 [TSl and rigorous mathematics 
|[T6l [TTl [T8l . While difficult to exactly achieve, various ap- 
proximate schemes make it more practical ll9l [T5l[T9ll20l . 

Exterior cloaking methods include cloaking due to anoma- 
lous resonance ETll22ll23ll24l where polarizable dipoles, and 
polarizable line quadrupoles, and clusters of arbitrarily many 
polarizable line dipoles in the vicinity of a flat or cylindrical 
superlens 11251 l26l 127 1 are cloaked, but apparently not larger 
objects [28 1, and cloaking due to complementary media l29l 
where an "antiobject" is embedded in a superlens, to create 
cancellation. 

Miller [ 10 1 found that active controls rather than passive 
materials could be used to achieve interior cloaking. Here we 
use active cloaking devices to achieve exterior cloaking. In 
principle this could be done by mimicking the effect of the 
cloaking device of ll29l . but our objective is to have an active 
cloaking device which does not require one to know the shape 
of the object to be cloaked. 

Perhaps an analogy with water waves can be made [14|. 
Our objective is to use the cloaking devices to create an area 
of still water, near but outside the cloaking devices, without 
disturbing the pattern of waves a certain distance away. Then 
a boat can be placed in the area of still water, without dis- 
turbing the surrounding waves: the boat is cloaked, and so 
are the cloaking devices. The area of still water is created by 
destructive interference between the surrounding waves and 
anomalously localized waves created by the cloaking devices. 



Our cloaking approach has the disadvantage that one needs to 
know in advance the incoming probing waves. 

For simplicity our analysis is restricted to the two dimen- 
sional case, corresponding to transverse electric or magnetic 
waves, so the governing equation is the Helmholtz equation. 
To begin with we study the two-dimensional quasistatic prob- 
lem since the analysis can be carried further in that case and 
yields valuable insights. For the Helmholtz equation our re- 
sults are purely numerical but provide convincing evidence 
that broadband exterior cloaking is possible. 

For the two dimensional quasistatic problem we assume 
that the dielectric constant of the background media is con- 
stant, i.e., the voltage is harmonic. Also by B r (£) C R 2 
we will denote the disk with radius r centered at x = (£,0), 
and by / we denote the potential due to exterior sources, that 
would exist in the absence of both the cloaking device and ob- 
ject to be cloaked. The exterior cloak we propose for the qua- 
sistatic problem, consists of one active cloaking device, which 
is a simply connected region containing the origin, along the 
boundary of which the potential can be prescribed, and a re- 
gion to be cloaked B a (S) with S > a > 0, that is exterior to 
the cloaking device. The cloaking is achieved as follows: ac- 
cording to the assumed known potential / (which could be ob- 
tained by suitably placing probes in the surrounding medium) 
the active device generates appropriate fields such that one 
will have very small fields inside B a (S) and at the same time 
a very small total scattering effect, (due to the device and 
the object), outside a sufficiently large disk, i.e., B 7 (0), with 
7 > a + S. Therefore, an exterior cloak, will, regardless of 
the probing field, achieve approximate invisibility of both the 
active device and any passive object placed in the disk B a (S). 

The question is then to provide a constructive way to find 
the potential at the device in order to achieve cloaking. First, 
notice that by subtracting / from the potential everywhere and 
applying an inversion transformation to the problem, namely 
z = 1/s where s — x\ + 1X2, the question now is to 

Find a harmonic v : R 2 — > R such that 

v « 0inB 1/7 (0) andv « — / in B at (5*). (1) 

Here we have a* = a/(S 2 ~ a 2 ) and 5* = 5/(8 2 — a 2 ). 
In fact we only need v to be harmonic outside the image of 



the cloaking device, but requiring it to be harmonic in all M. 2 
simplifies the problem. If such a function v exists then its 
Dirichlet data on the boundary of the cloaking device gives us 
the necessary potential one needs to generate at the surface of 
the device in order to achieve approximate invisibility. 

By introducing harmonic conjugate potentials one obtains 
the analytic extensions, V and F, of v and / respectively. 
Then problem ([TJ is equivalent to, 

Find V : C — > C analytic, such that 

V «0inB 1/7 (0) and V « -F in (2) 

Since the product of two analytic functions is again ana- 
lytic, the problem |2} can be equivalently formulated as 

Find W : C — > C analytic, such that 

IV « in B 1/7 (0) and ^ « 1 in B . (J t ). (3) 

To recover V one needs to multiply W by a polynomial which 
approximates — F in B Qt (5*). 

Next we consider the Hermite interpolation polynomial 
h : C — > C of degree 2n — 1 defined by 

h(0)=l, h(S,) = 0, 

ft(j')(Q) = fcO')(«y = q fori = l,...,n- 1. (4) 

From |4]), by algebraic and combinatoric manipulations to- 
gether with an induction argument, it can be shown that 




(5) 

Notice that |5]) implies the symmetry property h(8* — z) + 
h(z) = 1. 

Our claim that W(z) = 1 — h(z) satisfies the properties in 
<J3^> when I/7 and a* are small enough is strongly supported 
by fig. [T] the solid white curve corresponds to the contour 
\h{z)\ = O.Oland the dashed white curve to \h(z)-l\ = 0.01. 
Also, the last equality in (|5j implies that h(z) is in fact a power 
series in the variable , which by the ratio test converges as 
n — > 00 in the entire figure eight shaped region |z 2 — <5*z| < 
8 2 /4 (red curve in fig.[T|), and diverges everywhere outside this 
region, excluding the boundary. The convergence is uniform 
in any simply connected domain that lies strictly within the 
figure eight. Therefore the limit function is analytic in each 
half of the figure eight, and from the Taylor series of the limit 
function at the points z = and z = 8*, implied by Q, 
we deduce that the limit function is one in the left side of the 
figure eight, and zero in the right side of the figure eight. To 
ensure that <f3l) is satisfied we require that the disks _B 1 / 7 (0) 
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FIG. 1: Magnitude of the Hermite polynomial h(z) for <5„ = 1 and 
degree n = 10. The color scale is logarithmic from 0.01 (dark blue) 
to 100 (dark red). 

and B at (d*) lie, respectively, within the left and right sides of 
the figure eight, which is the case if both I/7 and a* are less 
than 5*/ (2 + 2\/2) . Notice that one gets cloaking (in the limit 
n — > 00) not just within the disk B ar (8*) but in the whole 
right half of the figure eight. Thus the scheme based on the 
Hermite polynomial defined in |5]) achieves exterior cloaking 
in the original s = Xjz plane for objects which are sufficiently 
small compared with 7. For example with <5* = 1, we have 
a* < 0.2, 7 > 4.9 and a = a*/(8 2 - a*) < 0.2. 

The cloaking device can be taken to be any simply con- 
nected region which contains the origin, but not the point 
(8,0). If it is a small disk, say B r (0) with r C 1 then one 
needs to generate enormously large potentials at its bound- 
ary. Alternatively its boundary could be taken as the con- 
tour where, say |/i(l/s)| = 100 (black curve in fig. [T] and 
El, then the cloaking device will tend to surround the cloak- 
ing region as n — > 00. Such a cloaking device is demon- 
strated in fig. [2] where we wish to cloak a (almost resonant) 
disk of radius 0.2 centered at s — 1.1 with dielectric constant 
e = —0.99 and located inside the cloaked region (solid white 
curve where |/i(l/s)| = 0.01). This "scatterer" deforms the 
vertical equipotential lines of the "incident" field f(s) = s 
when the device is inactive (fig. |2^). When the device is ac- 
tive the effect of the scatterer is greatly diminished, making it 
for all practical purposes invisible: outside the dashed white 
curve where \h(l/s) — 1| = 0.01, the equipotential lines are 
vertical (fig. |2j)). Indeed the discrepancy between the inci- 
dent field and the field on the dashed white circle in fig. |2j) is 
of about 1.1% of the incident field or 2.9% of the uncloaked 
scattered field, measured in the L 2 norm. Due to their associa- 
tion with high order multipoles these errors decay very rapidly 
as 7 increases. 

One possible extension of the quasistatic cloaking is to the 
Helmholtz equation Au + k 2 u = where k = 2tt/X is the 
wavenumber and A = 2irco /u> is the wavelength at frequency 
u> and constant propagation speed c . Numerical simulations 
suggest it is necessary to take D > 3 devices located at points 
Xi , X2 , . . . , X£> surrounding the cloaked region to design good 
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FIG. 2: Real part of the total field with the quasistatic cloaking device (a) inactive and (b) active in the s plane. The color scale is linear from 
— 10 (dark blue) to 10 (dark red) and the scatterer appears in black. 



cloaks regardless of the incident field direction. The cloaked 
region is for simplicity the disk |x| < a and the devices lie on 
the circle |x| = S. To ensure that objects inside the cloaked 
region are hard to observe at locations |x| > 7, the devices 
need to create a combined field Ud such that 

Ud ~ —Ui for |x| < a, and Ud ~ for |x| > 7, (6) 

where u t is the field due to exterior sources which would exist 
in the absence of the cloaking devices and the object to be 
cloaked. Since we want Ud to be a solution to the Helmholtz 
equation and decay far from the devices, we use the following 
ansatz, 

D N 

u d (x) = ^2 X! b m,nH^(k\^.-^. m \)exp[in9 m ], (7) 

m — 1 n——N 

where Hn is the n— th Hankel function of the first kind and 
9 m = arg(x — x m ) is the angle between the vectors x — x m 
and (1, 0). The coefficients 6 m „ 6 C are found numerically 
by enforcing ^ on N a points pf, . . . , of the circle |x| = 
a and 7V 7 points pj, . . . ,pj* of the circle |x| = 7. The 
resulting linear equations are Ab sa — and Bb sa 0, where 
b e C DM , M = 2N + 1, is a vector with the coefficients 
6 m , n , and the matrices A S C N <* xDM and B € C N -> xDM are 
constructed so that (Ab)j = u,i(p") and (Bb)j = Urf(pJ). 

Coefficients b satisfying these equations in the least squares 
sense can be obtained via the Singular Value Decomposition 
(SVD) in two steps. First a solution b with Ab ~ — U; 
is calculated using the truncated SVD. Then a correction z is 
found as a minimizer of ||B(bo + z)|| 2 such that Az = 0. The 
latter linear least squares problem can be easily solved with 
the truncated SVD if a basis of the nullspace of A is available, 
which is the case since the SVD of A was computed in the 
first step. Thus the coefficients to drive the cloaking devices 
are b = bp + z. 



We illustrate this procedure in fig. [3] with three devices lo- 
cated S = 10A away from the origin. Here k — c = 1, 
A = 2ir and we apply a plane wave incident field Ui(x) = 
exp[ikx ■ d] with d = (cos(27r/7), sin(27r/7)). The cloaked 
region is the solid white circle with radius a = 2A. Invisibility 
is enforced on the dashed white circle of radius 7 = 20A. The 
control points pj and p" are uniformly spaced and less than 
A/2 apart on their respective circles and N = 57 terms were 
used for the ansatz (|7]). Finally the scattered field resulting 
from the impenetrable "kite" obstacle with sound-soft (homo- 
geneous Dirichlet) boundary conditions is computed using the 
boundary integral equation approach in |30l . 

As shown in fig. when the cloaking devices are active 
they create a "quiet" region where the wave field is close to 
zero. An object lying in this region is practically invisible 
because both the scattered and device's fields are for all prac- 
tical purposes undetectable outside of the dashed white circle. 
The field outside the dashed white circle is nearly identical to 
the incident plane wave: the discrepancy on that circle is of 
about 1.5 x 1CT 5 % of the incident field or 6.8 x 1CT 5 % of 
the uncloaked scattered field, measured with the L 2 norm. Of 
course these relative errors depend on the particular choice of 
the cutoffs for the singular values in our two step approach. 

By Green's identities, these point devices could be re- 
placed by a bounded region containing xi , . . . , xc but not the 
cloaked region, provided the region's boundary has a control- 
lable single- and double-layer potential (i.e. a point source and 
dipole density 11301 ). One example is to take D disjoint disks, 
however if their radii are too small then the strength of the 
potentials on the disk boundaries would be enormous (since 
for n ^ 0, Hn\r) = C(r _ l"l) as r — > 0). A natural ques- 
tion is whether one can get exterior cloaking with reasonable 
field magnitudes. The black curves representing the contours 
|ud(x)| = 100 in fig. |3] suggest this is possible because they 
resemble three disjoint disks. 
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FIG. 3: Real part of the total field with the Helmholtz cloaking devices (a) inactive and (b) active. The color scale is linear from -1 (dark blue) 
to 1 (dark red) and the axis units are A. 



A complete mathematical discussion, including the rigor- 
ous analysis for the arguments for the quasistatic case, to- 
gether with a different analytical approach, will be included 
in a forthcoming publication. The theory for the Helmholtz 
problem remains an object of current research. We antici- 
pate that the results extend to three-dimensions and to the full 
Maxwell equations but this remains to be explored. 

The authors are grateful for support from the National Sci- 
ence Foundation through grant DMS-070978. 
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